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Пошук локальних екстремумтв у задач! пильного 
пакування неортентованих сфероконусв 


На основании квази Ф-функций построена математическая модель задачи плотной упаковки не- 
ориентированных сфероконусов в кубоиде минимальной высоты. Построенная модель позволила 
применить для поиска локальных экстремумов метод внутренней точки на последовательности под- 
областей области допустимых решений. Предложен метод построения разнообразных начальных точек. 
Представлен численный пример. 

Ключевые слова: квази Ф-функция, задача упаковки, сфероконус. 
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Кеу мога5: а4па$1 Ф-ВшсНоп, раскше рго ет, зрКегосопе. 


На основ! кваз1 Ф-функщй побудовано математичну модель задач! пильного пакування неортентованих 
сфероконус1в у кубо мимально! висоти. Побудована модель дозволила застосувати для пошуку локаль- 
них екстремум1в метод внутрииньс! точки на послдовност! шдобластей област! припустимих розв’язюв. 
Запропоновано метод побудови р1зноманиних початкових точок. Наведено числовий приклад. 
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Вопросы моделирования плотного размещения трехмерных объектов активно 
исследуются различными авторами во всем мире. Задачи, в которых приходится 
решать эти вопросы, широко используются в науке и промышленности [1-3]. Следует 
отметить, что на сегодняшний день разработано большое количество методов моде- 
лирования упаковок сфер. Однако в реальности очень часто приходится сталкивать- 
ся с задачами, в которых необходимо рассматривать трехмерные объекты произволь- 
ных пространственных форм. Моделирование плотного размещения таких объектов 
является ресурсоемкой задачей. В существующих работах по упаковке трехмерных 
объектов, форма которых отличается от сферической, рассматриваются только при- 
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ближенные подходы, не позволяющие провести адекватное моделирование задачи и 
получить хорошее решение. Один из распространенных подходов к решению задачи 
упаковки трехмерных объектов заключается в представлении объектов набором сфер и 
рассмотрении условий не пересечений этих сфер. Например, в работе [2] с помощью 
такого подхода решается задача упаковки цилиндров, сфероцилиндров, конусов, 
сфер и кубов. Еще один распространенный подход сводится к использованию 
простой «идеализированной» математической модели дискретизации объектов с 
помощью трехмерной сетки, ячейки которой являются кубоидами [3]. В работе [4] 
предлагается использовать сферополиэдры для представления различных трехмер- 
ных объектов. Используемый авторами алгоритм позволяет получать размещение 
большого числа объектов (>100) в задачах. Все эти подходы направлены на получе- 
ние быстрого, оценочного размещения большого количества объектов. 

Целью данного исследования является построение адекватной математи- 
ческой модели задачи плотного размещения сфероконусов и разработка метода 
поиска экстремумов данной задачи. 

Сфероконусом 5К называется выпуклый геометрический объект, который 
образован усеченным конусом К высоты 2 с радиусами верхнего и нижнего осно- 
ваний и и г; (н>1 ) соответственно, верхним сферическим сегментом С' с высотой и 


радиусом основания, равными у и и соответственно, и Нижним сферическим сегментом 
С, высота и радиус основания которого равны и и г соответственно. Условия, 


обеспечивающие выпуклость объекта имеют вид: 


ой и <, УВ) +4 (ив) 
› № -Р . 
(нь) +41? -("-")} 2 


Обозначим через и= (ирвмьмь,й) вектор метрических характеристик сферо- 


м и 


конуса. Сфероконус является «базовым» объектом, изменяя метрические характери- 
стики которого, можно получить следующие трехмерные объекты: сфероцилиндр 
(и= (Ирбьмь ий), п=Р); усеченный конус (и = (п,,,0,0,й); цилиндр (п = (1,^,0,0,), п=р); 
конус (и=(х,0,0,0,й))); сферические сегменты (и = (п,0,,0,0) или и = (0,0, %>,0)); 
сферический диск (= (ии, и, и», 0), п=Р). 

Пусть задано множество сфероконусов О; а — Г = а 

Сфероконусы О;1!=/, допускают аффинные преобразования трансляции на 
вектор и; = (х,,у;,2;)е ВЗ и поворота на углы 0; = (а,,В;)е В? вокруг осей координат 
Ох’ и О’у соответственно. Вектор движения объекта О;,Ёе[„ обозначим и; = (и;,6;). 


Тогда, вектор и=(щ,...,и„)е В“, ®, =5п, определяет положение всех объектов в 
трёхмерном пространстве. 
Пусть плоскость Г(Х,У) определяется уравнением Л(Х,У)=(М(У),Х)+4 =0, где 


У=(1,6,4)е В и М(У)= (51 х,—5т д с0$ 1,6036 с0$ 1). 


Плоскость г(х 2% ) делит пространство КЗ на два полупространства 


П- = (Хе? :Л(Х,У) < 0} и п* = КЕ ВЗ:л(х,У)>0]. 
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Область размещения Р (контейнер) описывается следующим образом: 


РЕК = (,у,2) КЗ: -а<х<а,-Ь<уЗЬь, н 525%} 


Контейнер Р высоты /й=й› -й обозначим Р(й). Тогда вектора у = (х›,6’, 4’) е К", 


ЕЕ 


= !1,...,6\, Порождающие плоскости, в которых лежат грани контейнера Р(й), имеют вид: 


О 


Задача. Найти вектор й Е №“, такой, что объекты О;, [е Г, , содержатся в кон- 
тейнере Р(й) и его высота й достигает минимума. 
Для построения математической модели поставленной задачи в работе ис- 
пользуется метод Ф-функций [5]. Эта математическая модель может быть пред- 
ставлена в виде: 


й* =штй, Х ЕЙ" < “+9742, (1) 
где 
ео, Ь]ЕГ1,1< }, 


ИИ =(ыу,н, ль) < К 69 (и 
Фи (и, 70) >0лЕ и Ей.2,.. ее. (2) 
Фу (муз, шах] (7, ,}, Ф,Р)= йа КЦ, 
Ибо ды быт ыы Я] 


вии.) Або, 7)+А, (в, М7) - а (в, м}. 


о2(и,,7)= ми а вом) 


А 


оз (и, Р)=^,,Р) *( А ад, умри, 


и’, 


2% 2% 


саб.) ло, Р-ьь (в, м пои (вез м®), 


рыб), | има Деб), 


752 — №52 


2. 2 2 2 
пон ед, 


2и.2 
В.(0,) = (т В, т и, соз В. ,с08 а, соз В}, 
— пп -1 
Ули ьаь у = + жьнбьь-ану „ву = ии 
Функция Ф, (... У®) является Ф-функцией объекта О;1е/[, и полупростран- 


ства П`, отделяемого плоскостью [х в , 

Использование метода Ф-функций позволило записать условия взаимного не- 
пересечения объектов в виде набора систем неравенств, левые части которых являются 
бесконечно дифференцируемыми функциями. 
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Задача (1), (2) является МР-трудной задачей. Для поиска приближения к гло- 
бальному экстремуму такой задачи необходимо разработать метод получения ло- 
кальных экстремумов. 

Поскольку математическая модель (1), (2) поставленной задачи построена в 
виде классической задачи нелинейного программирования, то для ее решения могут 
быть применены различные модификации методов нелинейной оптимизации (напри- 
мер, метод возможных направлений, метод внутренней точки и др.). 

Однако для применения численных методов нелинейной оптимизации не- 
обходимо иметь допустимую начальную точку. Среди методов, которые применяются 
для построения начальных точек, в задачах размещения объектов в основном исполь- 
зуются различные модификации «жадных» алгоритмов. Применение «жадных» алго- 
ритмов в рассматриваемой задаче представляет болыпую сложность. 

Поэтому, для получения начальных точек в работе предлагается метод, кото- 
рый предполагает, что все метрические характеристики объектов являются перемен- 
ными и принимают значения от некоторого минимального значения до изначально 
заданного. 

Рассмотрим предлагаемый метод решения. Для каждого объекта введём величину 


59, равную сумме его метрических характеристик. Таким образом, и = ме + 
ый и ы У + и ‚ Е[,. Также введём переменные величины „д; — коэффициенты гомо- 
тетии объектов О;, ГГ, и вектор к = (®,,...,К„)Е А”. 

Зададим высоту й = й° так, что объекты гарантированно размещаются в области 
Р(й°). Введём векторы и’ = (щ,...,и„)е В“, к, =(6,....=)е В", ки = < 
и сформируем точку Х’=(и’,к,,У’), где вектор и’ случайно генерируется таким 
образом, что и; ЕР] =, а вектор У’,, 
пендикулярную отрезку, соединяющему полюсы объектов О; и О ‚у БЕ < у, и 


.]Е[‚{</] ПОрождает плоскость, пер- 


проходящую через его середину. 
Сформулируем следующую задачу: 
п 
шах > ` 50, ХЕОс К^и тп, (3) 
1=1 
О = х = (им,к,У)с Ви +Фу+т ФА ри рАьК у, > 0,1, 7 Е Пл < У 


Ф, (к, 05 0,0 <; < Г, Ёе .2,.6 (4) 
Заметим, что существует такое = = (0/|, что Х"=(и',к,,У’)еО. Решим данную 


задачу, взяв точку Х’еО в качестве стартовой. Целью решения данной задачи 
является нахождение такого размещения объектов в области Р(п), при котором все 


метрические характеристики объектов достигают своих заданных значений, то есть 
К; =1, 1ЕГ„. В случае, если в результате решения задачи (3), (4) будет получена 


точка Х” [И кьУ й глобального максимума, то точка ("у "| берется в качестве 
начальной точки для поиска локального минимума задачи (1), (2). 

Для поиска локальных максимумов задачи (1), (2) и локальных минимумов 
задачи (3), (4) использовался метод последовательного улучшения значения функции 
цели на последовательности подобластей, образуемых наборами систем неравенств 
из квази Ф-функций для каждой пары сфероконусов. При этом поиск локального 
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экстремума на каждой подобласти осуществлялся с использованием библиотеки ГРОРТ [6]. 
В данной библиотеке реализованы алгоритмы внутренней точки, которые для поиска 
вектора движения используют матрицу Гессе левой части системы ограничений, 
описывающей область допустимых решений. 

Приведем численный пример. Положим количество объектов п=10, среди 
которых 4 «базовых» сфероконуса, 3 сфероцилиндра и 3 цилиндра. Размеры кон- 
тейнера установим а=25,Ь=15. Метрические характеристики и координаты точки 
локального минимума йеЙ’ приведены в табл. 1. Значение локального минимума: 


= 63,126. 


Таблица 1 — Метрические характеристики и координаты точки локального минимума 


о |0 
мм |] | Я [А | @ [| В 


ЗЕЕ пой 4996 
то 
0.785 | 1463 


1 

2 

3 

4 | Ь 

57 
6 

7 

8 

9 


.57Г [475 
2279 [2.555 
азот [ввв [988 о [о [131 [3.08 [ 103 [5.045 [0427 


10 | 3,251 | 6,281 | 6,281 0 0 21,57 | 8,493 | 19,94 | 5,25 | 4,709 


Рисунок 1 — Локальный минимум 


В работе благодаря использованию аппарата квази Ф-функций построена ма- 
тематическая модель задачи плотной упаковки сфероконусов в кубоиде. Построен- 
ная модель позволила применить для поиска локальных экстремумов метод внутрен- 
ней точки. 
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зеагстиз [ог Гос Емтета 
ор № п-Оптете4 эрйегосопеу Оепзе Расяия Ргоет 


ш {115 рарег а тафетайса| то4е| ог 4епзе раскше рго Мет оЁР поп-опещеа 
зрВегосопез шю а сибола 15 сопзнгисеа Бу изше Фе Ф-ЕЮпсНоп фестс. 

Ап аррПсаНоп ог Ф-ЕЮшсНопз$ аПо\уе4 из ю Юпишае тшиба| поп-пегзесноп$ 
сопа1юп$ Юг а рай оРоБ]ес{$ аз а зе оР тедиаНИез зузетз \шсь ей з14е5 аге шИпиеЙу 
ЧИегепнаЫе Рпсйоп$. О\лие ю 15 ЁсЕ а тафетайса! то4е! оЁ Ве ргоЫет 15 ргезеще4 аз 
а с1а551са| поп-Ппеаг ргоэгаплитте рго ет. 

Еог сопугасноп оЁ загИпе ро1ё5 \е ргорозе4 тео у\/№мсВ аззитез {На{ аП 
оБ]есё5 тейлс свагасег1$Ис$ аге уама е ап4 {1аКе уашез Нот а шшипиат опе {о Фе 
шаПу зресНед. 

п ог4ег 0 Впа 1оса| тшипа оР Фе ргоЫет (1), (2) ап4 оса] тахита ое рго ет 
(3), (4) ме изе Фе тефо4 оРсопипиочц$ пиргоуетеп: оЁ Пе об]есиуе ВапсНоп уае$ опа 
зедчепсе оЁ забгегтопз фа{ аге Югтеа Бу 5е{5 оР шедиаПиез зузбетз оРацаз1 Ф-ЕЮпсНоп$ 
Гог еасБ раш оЁ зрбегосопез. ТБе зеагсВ оЁа 10са| ехгетит ш еасВ забтезтоп 1$ рег- 
югтед Бу изте ТРОРТ Шгагу. 
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